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Abstract. Directional detection is a promising Dark Matter search strategy. Indeed, 
WIMP-induced recoils present a direction dependence toward the Cygnus constellation, while 
background-induced recoils exhibit an isotropic distribution in the galactic rest frame. Taking 
advantage on these characteristic features and even in the presence of a sizeable background, we 
show for the first time the possibility to constrain the WIMP properties, both from particle and 
galactic halo physics, leading to an identification of non-baryonic Dark Matter. However, such 
results need highly accurate track reconstruction which should be reachable by the MIMAC 
detector using a dedicated readout combined with a likelihood analysis of recoiling nuclei. 



1. Introduction 

Taking advantage of the astrophysical framework, directional detection of Dark Matter is an 
interesting strategy in order to distinguish WIMP events from background ones Indeed, hke 
most spiral galaxies, the Milky Way is supposed to be immersed in a halo of WIMPs which 
outweighs the luminous component by at least one order of magnitude. As the Solar System 
rotates around the galactic center through this Dark Matter halo, WIMPs should mainly come 
from the direction to which points the Sun velocity vector and which happens to be roughly in 
the direction of the Cygnus constellation (£q = 90°, Bq = 0°). Hence, we argue that a clear and 
unambigous signature of a Dark Matter detection could be done by showing the correlation of 
the measured signal with the direction of the solar motion. 



Several project of directional detectors are being developed and some of them are already 
taking data jH El |6] . The first step when analysing directional data should be to look for a signal 
pointing toward the Cygnus Constellation with a sufficiently high significance [7]. If no evidence 
in favor of a Galactic origin of the signal is deduced from the previous analysis, then an exclusion 
limit should be derived [8j. On figure [1] we present the projected discovery regions and exclusion 
limits for a forthcoming directional detector proposed by the MIMAC collaboration. We consider 
a 10 kg CF4 detector operated during ~ 3 years, allowing 3D track reconstruction, with a 10° 
angular resolution, a recoil energy range 5-50 keV and with a conservative background rate of 
10 evts/kg/year. From figure [H two different scenarii may be distinguished: 

• The dark and light-grey shaded areas represent the contours where a directional detection of 
Dark Matter would have a significance greater than 3a and 5a. These contours are deduced 




WIMP mass (GeV/c^ 



Figure 1. Spin dependent cross-section on proton (pb) as a function of the WIMP mass 
(GeV/c^). Exclusion limits from some direct detection experiments are presented, KIMS [9j and 
Picasso ^lOj as well as the theoretical region, obtained within the framework of the constrained 
minimal supersymmetric model from [11]. Contours corresponding to a significance greater 
than 3a and 5a are presented in dark and light grey. The exclusion limit corresponding to pure 
background data is presented as the black dashed line and the detector sensitivity as the black 
solid line. 



from the map-based likelihood method f7]. As an illustration, if the WIMP-nucleon cross- 
section is about 10~^ pb with a WIMP mass of 100 GeV.c"^, the detector would have a 
Dark Matter detection with a significance greater than 3a. 

• If the WIMP-nucleon cross-section is lower than 10^^ pb, then an exclusion limit is deduced 
using the extended likelihood method (black dashed line) presented in [8] . As a benchmark 
and to illustrate the effect of background on exclusion limits, the detector sensitivity (no 
event) is presented (black solid line). 

Figure [1] also presents exclusion limits from direct detection experiments, KIMS [9] and Pi- 
casso [To] as well as the theoretical region, obtained within the framework of the constrained 
minimal supersymmetric model taken from [11]. We can conclude that a directional detector 
like MIMAC will cover an important region of interest worth being investigated. 

Constraining Dark Matter parameters both from particle physics (mass rriy. and cross section 
an) and galactic halo (WIMP velocity distribution) with upcoming Dark Matter experiments is 
a main concern of current phenomenological studies. In the context of upcoming experiments 
which might give WIMP positive detection ((j„ ~ 10"'^ pb) instead of background rejection, it 
is of particular interest to develop a model-independent formalism to constrain Dark Matter 
properties. A new approach has been presented in [12]. The idea is to constrain Dark Matter 
properties using a high dimensional multivariate analysis within the framework of a general halo 
model. This strategy is referred to as quasi model-independent method as all parameters are 
constrained directly from the data of a single experiment. 

Indeed, directional detection presents a high identification potential thanks to the combined 
measurement of the recoil energy Eji and the recoil direction {iR,bfi). It allows to achieve 



multivariate analysis in order to obviate systematic bias in the determination of the WIMP 
properties, due to wrong halo model assumptions [H] . 



2. Identification of Dark Matter 

Within the framework of a quasi model-independent method and by using a Markov Chain 
Monte Carlo analysis of recoil events, we show for the first time the possibility to contrain the 
WIMP properties, both from particle physics {m^,an) and galactic Dark Matter halo physics 
(velocity dispersions). This leads to an identification of non-baryonic Dark Matter, which could 
be reached within few years by upcoming directional detectors p]. 



2.1. Directional framework 

Directional detection depends crucially on the WIMP velocity distribution. The isothermal 
sphere halo model is often considered but it is worth going beyond this standard paradigm in 
the case of a model-independent analysis. Indeed, recent results from N-body simulations are in 
favor of triaxial Dark Matter halos with anisotropic velocity distributions }13J. Moreover, recent 
observations of Sagittarius stellar tidal stream have shown evidence for a triaxial Milky Way 
Dark Matter halo [15] . 

The multivariate Gaussian WIMP velocity distribution [16j corresponds to the generalization 
of the standard isothermal sphere with a density profile p{r) oc l/r^, leading to a smooth WIMP 
velocity distribution, a flat rotation curve and no substructure. The WIMP velocity distribution 
in the laboratory frame is given by. 
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where a-y = diag[ax,(Ty,az] is the velocity dispersion tensor assumed to be diagonal in the 
Galactic rest frame (x, y, z) and Vq is the Sun motion with respect to the Galactic rest frame. 
When neglecting the Sun peculiar velocity and the Earth orbital velocity about the Sun, vq 
corresponds to the detector velocity in the Galactic rest frame and is taken to be = 220 
km.s~^ along the y axis pointing toward the Cygnus constellation at = 90°, 6© = 0°). The 
velocity anisotropy /3(r), is then defined [I7j as 
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According to N-body simulations, the (3 parameter at the Solar radius spans the range — 0.4, 
corresponding to radial anistropy. 

In the following, the input halo model used to generate simulated data, is chosen according 
to two models : a standard isotropic halo (/3 = 0) in which case the velocity dispersions are 
related to the local circular velocity vq = 220 km/s as cTx = cfy = o'z = ^^o/\/2; and an extremely 
anisotropic halo (/3 = 0.4), with the following velocity dispersions {ax = 200 km/s; cj^ = 169 
km/s; Gy = 140 km/s}. The latter case corresponds to the logarithmic ellipsoidal halo model 
from [16] with the Sun located on the major axis of the halo with the axis ratios p and q equal 
to 0.9 and 0.8 respectively. 



The directional recoil rate is given by [18] 
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with the WIMP mass, rrir the WIMP-nucleus reduced mass, po = 0.3 GeV/c^/cm^ the 
local Dark Matter density, ctq the WIMP-nucleus elastic scattering cross section, F(Eji) the 
form factor (using the axial expression from [19]) and fmin the minimal WIMP velocity required 
to produce a nuclear recoil of energy Er. Finally, f{vmm,Q) is the three-dimensional Radon 
transform of the WIMP velocity distribution f{v). Using the Fourier slice theorem [18j . the 
Radon transform of the multivariate Gaussian is, 
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In the following, the model is characterized by 8 free parameters which are 
{m-^,\ogiQ{an)-,lQ-,hQ,ax-,(yy,az-,Rb}, where the direction {lQ,bQ) refers to the main direction 
of the recorded events [7], (T„ is the WIMP-nucleon cross section directly related to fio in the 
pure proton approximation and Rf, is the background rate. We have considered flat prior for each 
parameter. In such case, the Bayes' theorem is simplified and the target distribution reduces to 
the 8 dimensional likelihood function ^{6) dedicated to unbinned data as, 
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where fig and pb = RbX(, are the expected number of WIMP and background events respectively. 
Rn refers to the energy and direction of each event while the functions S and B are the directional 
rate of the WIMP and the background events respectively. This MCMC analysis is based on the 
Metropolis-Hastings algorithm, using chain subsampling according to the burn-in and correlation 
lengths to deal only with independent samples [12j . 



2.2. Results from a benchmark model 

We first applied this MCMC method on a benchmark model with the following characteristics : 
isotropic halo model {(3 = 0) with a sizeable background contamination (10 kg^^.year"^), a 
50 GeV/c^ WIMP and a WIMP-nucleon axial cross section (j„ = 10~^ pb, well below current 
exclusion limits, see figure [TJ The input model is used to generate simulated data considering 
a MIMAC like detector as described above. These data are then analysed with the directional 
MCMC method. The full MCMC result with 2D correlations may be found in [12]. However, 
major results are shown on table [1] and the conclusion is threefold: 

• the discovery proof is given by the reconstruction of the main incoming direction (£0,60) 
pointing toward the Cygnus constellation within 2.5°, which is in favor of a positive Dark 
Matter detection. 

• The three velocity dispersions are strongly and consistently constrained according the input 
values corresponding to an isothermal sphere. Moreover, using equation ([2]), we can evaluate 
the posterior Probability Density Function (PDF) of the velocity anisotropy parameter (/3) 
which is in favor of an isotropic halo {(3 = — 0.073lJ^Q^g). 

• As the galactic Dark Matter halo properties are well constrained, constraints on the WIMP 
properties {rriy, an) are also strong and consitent according to the input WIMP parameters: 
(m^ = 5l.8tli^ GeV/c2, logio(cT„) = -3.0ll0;°| pb). 

As a conclusion of this benchmark study, we can deduce that this MCMC analysis, combined 
with directional data, is very well suited to constrain Dark Matter properties in a quasi model- 
independent way. 



Table 1. Comparison of the values of the parameters for the input model (isotropic halo and a 
flat background) and as extracted after the MCMC analysis from the marginalized distributions. 
We quote mean value of the PDF distribution and (68 % CL) error bars. 
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Figure 2. Left panel : 68% and 95% contour level in the (m-^,(T„) plan, for the isotropic input 
model and for a WIMP mass equal to 20, 50 and 100 GeV/c^. Right panel : posterior PDF 
distribution of the /3 parameter for the same models. 
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Figure 3. Left panel : 68% and 95% contour level in the (m-^, (T„) plan, for a 50 GeV/c^ WIMP 
and for two input models : isotropic (/3 = 0) and triaxial (/3 = 0.4). Right panel : posterior 
PDF distribution of the /3 parameter for the same models. 



2. 3. Varying the input parameters 

From the previous section, we have conclude that this MCMC analysis tool is very efficient to 
constrain Dark Matter properties. However, in this section we will briefly discuss the effect of 



some of the input parameters on the different constraints in order to estimate the performance 
of this analysis tool. We will first focus on the impact of the input WIMP mass. To do so, we 
have simulated three different sets of directional data corresponding to an input WIMP mass 
of = 20, 50, 100 GeV/c^ with a constant WIMP-nucleon cross-section an = 10~^ pb and the 
standard isotropic halo model. The results from the three MCMC runs are illustrated on figure 
El We present for the three WIMP masses, on the left panel, the 68% and 95% CL contours in 
the (m^^, log]^o((Jn)) plan and on the right panel, the posterior PDF P{/3\D) of the anisotropy 
velocity parameter (3. 

It can be deduced from figure [2] that the constraints strongly depend on the input WIMP mass, 
but in each case, the constraints are consistent with the input values. Then, this analysis has 
been shown to be working for any input WIMP mass although the constraints are stronger for 
light WIMPs. This is due to the fact that the signal characteristics, i.e the slope of the energy 
distribution and the width of the angular distribution, evolve slowly with the WIMP mass once 
> 100 GeV/c^, as shown in ff\. 

In the following, we investigate the effect of an extremely triaxial halo model with f3 = 0.4 
on the estimation of the Dark Matter parameters (my., an,/3)- The results from the MCMC run 
on a simulated dataset corresponding to a WIMP mass of 50 GeV/c^ with the anisotropic halo 
model are presented on figure [3l As previously, on the left panel is presented the constraint at 
68% and 95% on the {m^,logiQ{an)) plan while on the right panel is given the deduced posterior 
PDF of the f3 parameter. For convenience and comparison, the results from the benchmark 
input model (isothermal sphere with a 50 GeV/c^ WIMP) are recalled. 

Prom the left panel of figure [3l we can conclude that the two halo models give similar constraints 
which are both consistent with the input values. In fact, and as foreseen, the fact that the ve- 
locity dispersions are set as free parameters in the MCMC analysis allows to avoid induced bias 
due to wrong model assumption. From the right panel of figure [3] we can deduce that the /3 
parameter is well constrained: /3 = 0.381q'^ and strongly in favor of an anisotropic Dark Matter 
halo. 

As a conclusion of this study, it should be highlighted that the combination of information 
from the angular and energy distributions leads to robust allowed regions in the {m^,logiQ{an)) 
plan, since the halo model is also being constrained with the MCMC analysis from the same 
dataset of a single directional detection experiment. Moreover, the velocity anisotropy parameter 
/3, i.e. the three velocity dispersions, could be sufficiently constrained to discriminate between 
different halo models with future directional detectors such as the one proposed by the MIMAC 
collaboration j^. 

3. Track reconstruction vi^ith MIMAC 
3.1. MIMAC prototype 

The MIMAC prototype is the elementary chamber of the future large matrix. It allows the 
possibility to show the ionization and track measurement performance needed to achieve the 
directional detection strategy. The primary electron-ion pairs produced by a nuclear recoil in 
one chamber of the matrix are detected by drifting the electrons to the grid of a bulk micromegas 
[20] and producing the avalanche in a very thin gap (128 or 256/im). The electrons move towards 
the grid in the drift space and are projected on the anode thus allowing to get information on 
X and Y coordinates. To access the X and Y dimensions with a 100 fim spatial resolution, a 
bulk micromegas with a 4 by 4 cm^ active area, segmented in pixels with a pitch of 350 //m is 
used as 2D readout. In order to reconstruct the third dimension Z of the recoil, a self-triggered 
electronics has been developed. It allows to perform the anode sampling at a frequency of 40 
MHz. This includes a dedicated 16 channels ASIC [21] associated to a DAQ [22]. The total 



Table 2. Comparison between the input and the output values of the track characteristics with 
the {S = downward} hypothesis for a 100 keV proton recoil in 50 mbar of pure isobutane. We 
quote mean value of the PDF distribution and (68 % CL) error bars. 



X [cm] Y [cm] Z [cm] 9 [°] (p [°] 

Input 5 45 

Output -0.003 ±0.009 0.006 ±0.012 5.068 ±0.18 45.47 ±1.7 -0.59 ± 2.6 



recoil energy is deduced from the measured ionization quenching factor (IQF) [23]. 

In order to fully exploit the data from the MIMAC detector and to constrain Dark Matter 
properties, the recoiling tracks have to be accurately reconstructed: direction, sense and position 
in the detector volume. To achieve this goal, we first developed a track simulation software using 
SRIM [23] (simulated tracks), Magboltz [25] (electron drift velocity and dispersions) combined 
with the detector response simulation (MIMAC DAQ). This way, we can simulate any kind 
of recoiling track according to the following input parameters: X, Y, Z, 9, (f) and the sense S 
(upward or downward). Each track is characterized by three types of topological observables: 
number of slices, position of the center of gravity X^°'^y and y''"'^'^ of the collected charge for 
each slice and the width AX^°'^y and AY^"'^'^ along X and Y of each slice. Hence, the number 
of topological observables associated to any recoil track is given by: 

Nobs = l+^X N slice (6) 

In this paper, we present a new MIMAC data analysis, based on [26j . in order to recover the 
track parameters: {X, Y, Z, 0, (/>, 5}. As a working example, we will consider in the following a 
simulated track of an hydrogen recoil of 100 keV at the position {X = 0, Y = 0, Z = 5} cm, in 
the direction {9 = 45°, = 0°} and going downward. The analysis is based on the computation 
of the likelihood function ^{X, Y, Z, 9, (j)\S) defined as, 

^{X, Y, Z, 9, ^\S) = P{NsUce) n ^(^^'')^'(i"^''')^(AX^^^)P(Ay„^-^) (7) 

71=1 

where the different probability P are estimated using a frequentist approach by simulating a 
large number of tracks for each step in the parameter space. 

In order to compute the likelihood function, we used a Markov Chain Monte Carlo algorithm 
and the full MCMC result following the {S = downward} hypothesis is given in [26]. Deduced 
constraints on the five parameters and comparison with their input values are described in ta- 
ble [2j From this single track analysis it can be deduced that the spatial uncertainty is of the 
order of 100 fim for X and Y and about 2 mm for the absolute height Z; and the angular un- 
certainty is about 3°. However, even if systematic studies are needed to accurately estimate the 
different resolutions, results from table[2]are strongly in favor of accurately reconstructed tracks. 

So far, we did not discuss the possibility to discriminate between the two hypothesis: 
downward or upward, the so-called "head-tail effect". To do so with our track reconstruction 
method, we compute the likelihood function according to the second hypothesis {S = upward}. 
Then, a direct comparison of the maximum Log-Likelihood values of each hypothesis (-3.8 for 



{S = downward} and -5.1 for {S = upward}) allows to conclude in favor of a downwarding 
track consistent with the input simulated track. Hence, this method seems to be able to achieve 
the sense recognition of the recoiling tracks, but here again, systematic studies must be done to 
estimate its efficiency [26]. 

4. Conclusion 

In this paper, we have shown that the directional detector proposed by the MIMAC collaboration 
should be sensitive to an important region of interest of the Dark Matter parameter space {m^, 
an) motivated by constrained minimal supersymmetric models. In the case of a strong WIMP 
positive detection with (T„ ~ 10~^ pb, an identification of Dark Matter could be achieve in a 
model independent way using the MCMC analysis method presented in section 2. We have shown 
that constraining the velocity dispersions of the galactic Dark Matter halo allows to obviate bias 
in the determination of the WIMP properties, due to wrong halo model assumptions. Moreover, 
we have shown the possibility to discriminate between various halo model characterized by 
their velocity anisotropy parameter. However, such interesting results requires accurate track 
reconstruction: direction, position and sense. To reach this goal, we developed a new track 
reconstruction method based on a likelihood analysis which has shown very promising results, 
i.e high resolutions. Although, systematic studies must be done to accurately estimate the 
differents resolutions and the sense recognition efficiency. 
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